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I.  Introduction 


Starting  in  September  1992,  the  DoD  AASERT  program  has  provided  a  three  year 
augmentation  award  (AFOSR-F49620-92-J-0349)  to  the  PI  for  training  two  graduate 
students  to  perform  scientific  research  in  conjunction  with  the  Pi's  ongoing  research 
activities:  "Scattering,  Absorption  and  Generation  of  Microwaves  in  time  varying  and 
structured  plasmas"  then  sponsored  by  grant  AFOSR-9 1-0002,  and  currently  by  the  grant 
AFOSR-F49  670-94-0076.  During  the  past  three  years,  the  two  students  Mr.  James  Faith 
and  Ms.  Ruthie  Lyle  have  made  significant  progress  in  both  academic  area  (coursework) 
and  research  area  (carrying  out  research  projects)  for  their  thesis  dissertations.  Both  have 
passed  the  area  examination,  and  have  completed  all  required  coursework,  are  now  in  the 
final  stage  of  their  Ph.D.  thesis  work.  Contained  in  this  report  are  the  works  completed  so 
far.  Their  theses  representing  the  completed  final  report  will  be  finished  in  a  few  months 
and  will  be  submitted  to  amend  the  present  report. 

The  two  students'  research  involves  two  aspects  of  microwave-plasma  interaction. 
The  area  of  investigation  of  Mr.  James  Faith  is  primarily  on  with  the  scattering  and 
generation  of  microwaves  by  a  spatially  periodic  and  rapidly  time  varying  plasma.  The 
plasma  is  produced  by  a  nanosecond  (nS)  dc  discharge  through  six  pairs  of  equally  spaced 
parallel  electrodes.  The  results  indicate  that  this  plasma  device  works  as  a  dc  to  ac 
converter.  A  wideband  microwave  pulse  is  produced  during  each  discharge.  The  main 
goal  of  this  work  is  to  use  such  a  fast  growing  plasma  to  upshift  the  frequency  of  an 
incident  CW  microwave  from  an  external  source.  It  turns  out  that  an  unexpected 
downshift  process  is  observed.  The  downshifted  line  is  trapped  in  the  periodic  plasma  in  a 
way  similar  to  the  trapping  of  an  optical  signal  in  a  band  gap  of  a  periodic  dielectric 
medium.  Computer  simulations  have  confirmed  this  observation.  James  has  also  worked 
on  nonlinear  wave  particle  interaction  which  leads  to  chaos  in  the  particle  trajectory.  This 
work  is  enclosed  in  the  appendix. 

Ms.  Ruthie  Lyle  (a  very  motivated  African  American  woman,  a  very  rare  candidate 
in  the  engineering  area)  is  the  other  student  supported  by  the  DoD  AASERT  program. 
Her  research  deals  with  the  potential  application  of  plasma  electromagnetic  absorption  and 
scattering  for  shielding  objects  from  detection  by  conventional  radar  systems.  A  quasi¬ 
particle  theory  based  on  Wigner  distribution  function  has  been  developed  to  study  wave 
propagation  and  scattering  in  a  periodically  structured  plasma  This  theory  includes  the 
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effect  of  multiple  scattering  and  is  applicable  in  the  cases  of  interest.  Ruthie  has  applied 

the  theory  to  study  the  ionospheric  scintillation  phenomena.  Her  analysis  includes  the 

spectral  bandwidth  effect  of  ionospheric  density  irregularities  on  the  scintillation  results. 

II.  Publications 

Journal  Articles: 

1.  James  Faith,  S.P.  Kuo,  and  Joe  Huang,  "Chaotic  Electron  Motion  Driven  by  Whistler 
Waves  in  the  Magnetosphere",  accepted  for  publication  in  Comments  on  Plasma 
Physics  and  Controlled  Fusion. 

2.  James  Faith,  S.P.  Kuo,  and  Joe  Huang,  "Electron  Precipitation  Caused  by  Chaotic 
Motion  in  the  Magnetosphere  due  to  Large  Amplitude  Whistler  Waves",  submitted 
to  Journal  of  Geophysical  Research. 

Proceedings  Issued  Articles: 

1.  James  Faith,  S.P.  Kuo,  and  Joe  Huang,  "Chaotic  Electron  Motion  Driven  by  Whistler 
Waves  in  the  Magnetosphere",  International  Conference  on  Phenomena  in  Ionized 
Gas,  August,  1995,  Hoboken,  NJ. 

2.  Ruthie  Lyle,  "The  Spectral  Effect  of  the  Ionospheric  Irregularities  on  the  Scintillation 
of  Transionospheric  Signals",  International  Conference  on  Phenomena  in  Ionized 
Gas,  August,  1995,  Hoboken,  NJ. 

Conference  Presentations: 

1.  J.  Kim,  J.  Faith,  S.P.  Kuo,  "Design  of  a  High  Harmonic  Cusptron  Oscillator  of  an 
Inverted  Configuration" ,  APS  Annual  Meeting  of  Plasma  Science,  11/16-11/18,  1992, 
Seattle,  WA. 

2.  J.  Kim,  J.  Faith,  and  S.  P.  Kuo,  "Numerical  Simulation  for  High  Harmonic  Inverted 
Cusptron  Devices",  IEEE  International  Conference  on  Plasma  Science,  611-619, 
1993,  Vancouver,  Canada. 

3.  J.  Faith,  J.  Kim,  S.P.  Kuo,  "Design  of  High  Harmonic  Inverted  Cusptron  Device", 
APS  Annual  Meeting  of  Plasma  Science,  11/1-11/5,  1993,  St.  Louis,  MO. 
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4.  S.P.  Kuo,  J.  Kim,  J.  Faith,  "Design  of  High  Harmonic  Inverted  Cusptron  Device", 
DoD  Advisory  Group  on  Electron  Devices  Conference,  5/10-5/12,  1994,  Monterey, 
CA. 

5.  J.  Faith  and  S.P.  Kuo,  "Generation  of  ELF  and  VLF  Waves  in  the  HF  Heater- 
Modulated  Polar  Electrojet",  IEEE  International  Conference  on  Plasma  Science,  616- 
6/8,  1994,  Sante  Fe,  NM. 

6.  J.  Kim,  J.  Faith,  S.P.  Kuo,  "High  Harmonic  Cusptron  Device  with  an  Inverted 
Structure",  IEEE  International  Conference  on  Plasma  Science,  6/6-6/8,  1994,  Sante 
Fe,  NM. 

7.  J.  Faith,  S.P.  Kuo,  J.  Sun,  J.  Huang,  "Experimental  Study  of  Frequency  Shift  of  an 
Electromagnetic  Wave  Propagating  Through  a  Rapidly  Created  Periodic  Plasma",  APS 
Annual  Meeting  of  Plasma  Science,  1 1/7-1 1/1 1, 1994,  Minneapolis,  MN. 

8.  J.  Sun,  J.  Huang,  S.P.  Kuo,  J.  Faith,  "Simulation  Study  on  Electromagnetic  Frequency 
Shift  in  Rapidly  Created  and  Periodically  Structured  Plasma",  APS  Annual  Meeting  of 
Plasma  Science,,  11/7-11/11,  1994,  Minneapolis,  MN. 

9.  Ruthie  Lyle,  J.  Huang,  and  S.P.  Kuo,  "Electromagnetic  Wave  Scattering  in  an 
Inhomogeneous  Plasma  Characterized  by  Periodic  Structure",  APS  Annual  Meeting  of 
Plasma  Science,,  1 1/7-1 1/11, 1994,  Minneapolis,  MN. 

10.  J.  Faith,  J.  Huang,  S.P.  Kuo,  "Experimental  Study  of  Frequency  Shift  of  an 
Electromagnetic  Wave  Propagating  Through  a  Rapidly  Created  Periodic  Plasma  , 
IEEE  International  Conference  on  Plasma  Science,  6/5-6/S,  1995,  Madison,  WI. 

11.  J.  Faith,  J.  Huang,  S.P.  Kuo,  "Experimental  Study  of  the  Conversion  of  DC  Electric 
Fields  to  Microwave  radiation  by  an  Ionization  Front  Created  by  Successive 
Discharges"  IEEE  International  Conference  on  Plasma  Science,  6/5-6/S,  1995, 
Madison,  WI. 

12.  J.  Huang,  J.  Faith,  S.P.  Kuo,  "Chaotic  Particle  Motion  in  Large  Amplitude  Whistler 
Wave  in  the  Magnetosphere",  IEEE  International  Conference  on  Plasma  Science, 
6/5-6/S,  1995,  Madison,  WI. 

13.  Ruthie  Lyle,  S.P.  Kuo,  and  J.  Huang,  "Study  of  Electromagnetic  Wave  Scattering  by 
Periodic  Density  Irregularities  in  Plasma",  IEEE  International  Conference  on  Plasma 
Science,  6/5-6/S,  1995,  Madison,  WI. 

14.  J.  Huang,  S.P.  Kuo,  and  James  Faith,  "Multiple  Frequency  Conversion  of  an 
Electromagnetic  Wave  Interacting  with  a  Rapidly  Created  Periodic  Plasma",  APS 
Annual  Meeting  on  Plasma  Science,  11/6-11/10,  1995,  Louisville,  KY. 
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15.  S.P.  Kuo,  James  Faith,  and  J.  Huang,  "Frequency  Shift  of  an  Electromagnetic  Wave 
Interacting  with  a  Rapidly  Created  Periodic  Plasma",  APS  Annual  Meeting  on  Plasma 
Science,  11/6-11/10, 1995,  Louisville,  KY. 

16.  J.  Faith,  Joe  Huang,  and  S.P.  Kuo,  "Experiment  on  the  Conversion  of  dc  Electric 
Fields  to  Radiation  by  Successive  Discharges",  APS  Annual  Meeting  on  Plasma 
Science,  11/6-11/10,  1995,  Louisville,  KY. 

17.  James  Faith,  S.P.  Kuo,  Joe  Huang,  "Chaotic  Electron  Motion  Driven  by  Whistler 
Waves  in  the  Magnetosphere",  APS  Annual  Meeting  on  Plasma  Science,  1 1/6-1 1/10, 
1995,  Louisville,  KY. 

18.  Ruthie  Lyle,  S.P.  Kuo,  and  J.  Huang,  "Spectral  Effect  of  Ionospheric  Irregularities  on 
Scintillation  of  Transionospheric  Signals",  APS  Annual  Meeting  on  Plasma  Science, 
11/6-11/10,  1995,  Louisville,  KY. 

19.  James  Faith,  S.P.  Kuo,  and  Joe  Huang  /'Electron  Precipitation  caused  by  Chaotic 
Electron  Motion  in  the  Magnetosphere  due  to  Large  Amplitude  Whistler  Waves", 
National  Radio  Science  Meeting,  1/9-1/13,  1995,  Boulder,  CO. 

20.  Ruthie  Lyle,  S.P.  Kuo,  and  J.  Huang,  "Spectral  Effect  of  Ionospheric  Irregularities  on 
the  Scintillation  of  Transionospheric  Signals",  National  Radio  Science  Meeting,  1/9- 
1/13,  1995,  Boulder,  CO. 


HL  Summary  of  the  Works  Completed 

A.  Scattering  of  Microwaves  by  a  Rapidly  Produced  Spatially  Periodic  Plasma 

Recently,  there  has  been  a  considerable  interest  in  the  study  of  the  propagation  of 
an  electromagnetic  wave  through  a  plasma  whose  density  profile  changes  rapidly  in  time. 
It  is  motivated  by  the  fact  that  a  plasma  device  may  be  devised  to  tune  the  frequency  of 
existing  sources  of  microwaves  or  to  convert  a  CW  source  wave  into  a  wideband  pulse 
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train,  where  the  period  of  the  pulse  train  and  the  length  of  each  pulse  are  determined  by 
the  repetition  rate  of  the  plasma  discharge  and  the  length  of  the  plasma  volume 
respectively.  Such  a  wideband  pulse  train  has  a  potential  application  as  a  radar  jamming 
source 


Theory  and  numerical  simulation  on  this  wave  frequency  shift  by  interacting  with  a 
suddenly  created  uniform  plasma  have  been  performed.  Several  experiments  confirming 
the  theoretical  predictions  have  been  carried  out.  In  those  works,  a  plasma  slab  (with 
plasma  frequency  co^)  is  used  to  convert  the  preexisting  electromagnetic  wave  (with 

frequency  ©0)  to  an  upshifted  frequency  (©o  +©2C)1/2.  The  physical  mechanism  can  be 

simply  explained  as  follows.  Since  the  plasma  slab  is  generated  suddenly  in  time,  the  wave 
does  not  experience  any  spatial  variation  inside  the  slab,  and  therefore  the  wavelength  of 
the  wave  will  not  change.  The  wave  in  the  new  medium  has  to  satisfy  the  new  dispersion 

relation  ©2  =®2pe  +  klc2  and  k0  =  ©0/c,  thus  its  frequency  ©  is  upshifted  from  the  initial 
value  ©0  =  k0c  to  the  new  value  ©  =  (©j;  +©2C),/2 

In  the  present  work  we  consider  the  case  of  a  wave  propagating  through 
periodically  spaced  suddenly  generated  plasma  slabs.  The  frequency  shift  phenomena  in 
this  case  becomes  much  more  complicated  since  the  wave  k0  can  now  interact  with  the 

structure  ^of  the  plasma  and  different  Floquet  modes  of  the  periodic  structure  can  be 

excited.  Therefore,  a  number  of  frequencies  are  expected  to  emerge  due  to  the  temporal 
discontinuity  a  periodically  spatial  discontinuity  the  wave  experiences.  The  result  of  this 
scattering  phenomena  can  be  concisely  represented  in  a  band  diagram  of  the  dispersion 
curve  similar  to  that  encountered  in  solid  state  physics.  This  curve  compared  to  that  for  a 
uniform  plasma  shows  that  pass  bands  emerge  in  the  cutoff  frequency  regime  while  the 
stop  bands  (band  gaps)  appear  in  the  propagating  frequency  regime.  Therefore  both 
upshifted  and  downshifted  signals  are  expected  to  be  generated  as  the  initial  wave  interacts 
with  such  a  suddenly  created  periodic  plasma.  The  experiment  has  been  performed  to 
confirm  the  theoretical  predictions. 

We  have  constructed  a  chamber  experiment  to  investigate  the  frequency  shifting  of 
an  CW  microwave.  The  experiment  consists  of  a  pair  of  electrodes,  each  consisting  of  an 
aluminum  sheet  in  which  eleven  3cm  slots  have  been  cut,  with  3cm  between  each  slot. 
These  are  backed  with  a  solid  copper  sheet  which  provides  an  effective  waveguide  for  the 
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incident  wave.  In  each  of  the  depressions  created  by  the  slots  in  the  aluminum  we  arrange 
many  metal  spikes.  The  discharge  will  occur  only  in  the  region  where  the  spikes  are 
present  (as  opposes  to  arcing  from  the  comers  of  the  electrodes),  due  to  the  large  electric 
fields  at  the  tips.  A  photograph  of  the  discharge  is  shown  in  Fig.  1.  These  electrode 
stmctures  are  placed  facing  one  another  within  a  vacuum  chamber  (dimensions: 
2Vi'xl,A'xT)  constructed  of  1"  Plexiglas.  When  evacuated  to  low  pressure  (~0.5  Torr)  we 
rapidly  create  a  plasma  between  the  electrodes  by  discharge  from  a  Marx  capacitor  bank. 
The  CW  microwave  is  fed  into  one  end  of  the  chamber  with  an  S  band  horn  antenna,  and 
the  received  signal  is  collected  by  another  hom  of  the  appropriate  band  at  the  other  end. 
The  results  of  several  hundred  discharges  are  collected  on  our  HP  8105  spectrum  analyzer 
and  downloaded  to  a  PC  for  storage.  A  schematic  of  the  experimental  setup  is  shown  in 
Fig.  2. 


We  have  investigated  two  cases  experimentally.  In  the  first,  the  incident  frequency 
is  chosen  to  be  4GHz,  and  the  discharge  75  kV,  such  that  ©>©p.  In  this  case  the  dominant 
spectral  features  are  up-shifted,  to  be  discussed  in  the  theory  below.  We  find  an  upshifted 
peak  at  7.5  GHz,  which  is  a  greater  relative  shift  than  that  obtained  by  a  single  slab 
plasma.  This  result  is  presented  in  Fig.  3.  It  should  be  noted  that  in  this  case  the 
conversion  efficiency  is  very  low.  In  the  second  experimental  case  investigated,  incident 
frequency  of  3.05  GHz  and  a  100  kV  discharge  with  ©<©p,  we  encountered  some  new 
results,  unexpected  for  a  uniform  plasma  slab.  We  find  that  the  dominant  new  wave  can 
be  trapped  between  the  plasma  layers  and  its  downshifted  peak  at  2.5  GHz.  As  shown  in 
Fig.  4.  What  is  noted  that  there  is  a  tremendous  improvement  in  the  conversion  efficiency, 
with  the  downshifted  peak  1 5  dB  down  from  the  incident  wave.  This  represents  about  a 
30  dB  improvement  over  the  largest  upshifted  peak  we  have  observed  (e.g.  Fig.  3). 

The  experimental  results  can  be  understood  qualitatively  by  using  the  Floquet 
theory.  It  is  done  by  considering  a  wave  of  the  type: 

f/fexp[/£0(l-©2 /©2)1/2z]+2?exp[-/£0(l-©2 /©2)1/2z]  inside  the  slabs 

*  “  j  ) 

|C exp[/£0 (z-  L)]  +  D exp[-/£0 (z-L)]  outside  the  slabs 

where  L  is  the  periodicity,  d  is  the  slab  width,  and  k0  =  2n/X0  is  the  free  space 
wavenumber.  Applying  the  boundary  conditions  on  'F  and  its  derivative,  along  with  the 
Bloch  wave  condition: 

¥(z)  =  exp[-/pZ]'F(z  +  Z) 

we  can  find  the  dispersion  relation  for  the  wave  propagating  in  the  periodic  plasma  as: 
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Figure  1  A  photograph  of  the  discharge. 
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Figure  2  A  schematic  of  the  experimant. 
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Figure  3  The  experimantally  observed  down  shifted  peak. 


cos (PZ,)  =  cos(k0r\d)  cos [k0(d  -  L )] 

+(1  /  2)(r|  + 1  /  rj)  sin(^ri^)  sin[£0(</  - 1)]  (2), 

where  r|  =  (l-©2 /©2)1/2  the  effective  dielectric  constant  of  the  plasma  regions,  and  P  is 

the  propagation  constant  for  the  periodic  plasma.  Shown  in  Fig.  5  is  an  example  of  the 
band  diagram.  An  examination  of  the  regions  where  P  is  real  gives  the  regions  where  we 
would  expect  to  find  shifted  peaks.  Fig  5  indicates  that  a  number  of  small  pass  bands 
appear  in  the  region  ©  <©J),  consistent  with  the  experimental  observations  (Fig.  4). 

This  estimate  of  the  output  spectrum  can  be  confirmed  using  a  finite  difference 
time  domain  method  to  solve  the  wave  equation: 

\  <(*>')  E(z  _  0  (3) 
dz 2  c2  dt2  c 2  ’ 

where  ©p  is  a  function  of  time  with  a  periodic  spatial  profile.  After  solving  for  the  time 

evolution  of  the  field,  we  apply  a  Fast  Fourier  Transform  to  the  time  series  to  obtain  the 
spectrum.  These  simulations  confirm  the  predictions  of  the  dispersion  relation  and  show: 

1 .  when  dlL  is  small  many  Floquet  modes  can  be  excited,  with  many  upshifted  peaks, 
(Fig  6). 

2.  when  ©^  /  ©0  is  large  we  can  obtain  frequency  components  greater  than  (©o  +  ©2)1/2, 
or  downshifted  to  ©  <  ©  p  (Fig.  7). 

3.  when  d /  L  « 1  the  peak  frequencies  all  approach  one  another  near  (©2  +©2  )1/2 . 

Thus,  the  periodic  plasma  offers  many  new  interesting  phenomena  over  the  case  of  a  single 
slab. 
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Figure  5 

The  o) -p  dispersion  relation  (or  band  diagram)  L  =  Xq,  d/L  =  0.6,  G)p  /  coq  =  1.2 


to 


The  dispersion  curve  exhibit  a  band  structure  in  the  frequency  range  corresponding  to  an 
underdense  plasma.  While  in  the  frequency  range  corresponding  to  an  overdense  plasma  (below  the 
plasma  frequency),  small  pass  bands  emerge. 
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When  d  /  L  is  small,  many  Floquet  modes  can  be  excited.  Moreover,  if  cop  /  o)q  is  large,,  the 
periodic  plasma  can  effectively  shift  the  initial  wave  frequency  to  frequencies  larger  than 


Figure  7 

A  numerical  simulation  of  the  problem,  showing  similar  spectral 
features  to  those  observed  experimentally. 


B .  Study  of  Electromagnetic  Wave  Scattering  in  Periodically 

Structured  Plasma 

Radar  cloaking  requires  a  significant  attenuation  to  the  backscattered 
radar  signal.  There  are  two  approaches  to  achieve  this  requirement.  One 
approach  is  to  introduce  an  absorbing  material  in  front  of  the  target,  the  other 
to  introduce  a  structured  material  instead.  The  later  situation  provides  for 
significant  side  scattering  of  the  incident  signal  and  is  believed  to  be  more 
effective  and  feasible  for  the  radar  cloaking  purpose.  Therefore,  the  present 
work  is  devoted  to  the  development  of  a  theory  describing  wave  scattering 
phenomenon  in  a  periodically  structured  plasma. 

The  convention  scattering  theory  is  based  on  forward  scattering 
approximations  and  thus,  is  limited  to  small  density  perturbation.  As  the 
density  perturbation  increases,  multiple  scattering  becomes  important.  If  the 
perturbation  expansion  is  convergent,  the  conventional  scattering  theory  can 
still  be  applied,  but  a  tedious  iteration  procedure  is  involved.  Therefore,  a 
quasi-particle  theory  which  includes  the  multiple  scattering  process  in  its 
formulation  is  developed.  The  propagating  wave  is  treated  as  a  group  of 
quasi-particles  whose  density  and  momentum  distribution  are  proportional  to 
the  spatial  and  spectral  distribution  of  the  wave  intensity,  respectively. 
Hence,  the  distribution  of  quasi-particles  in  the  phase  space  can  be 
characterized  by  a  Wigner  distribution  function  which  is  governed  by  a 
"Boltzman"  type  equation  with  an  effective  collision  term  arising  from  the 
quasi-particles  interaction  with  plasma. 

This  theory  is  applied  for  explaining  our  laboratory  experiment,  as  well 
as  the  experimental  observations  in  space  plasmas  such  as:  the  scintillation  of 
Beacon  satellite  signals  and  ionospheric  sounding  (i.e.  spread-F).  The  work 
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presented  in  the  following  is  focused  on  the  study  of  scintillation  phenomena. 
The  scintillation  of  transionospheric  signals  is  an  area  of  current  interest  to 
the  Geophysical  directory  of  the  Air  Force  Phillips  Laboratory.  The 
conventional  theories  are  not  applicable  for  strong  scintillation  (i.e. 
Scintillation  index  S4  >  0.3).  However,  many  experimental  observations 
conclude  that  S4  >  0.6  is  not  unusual  in  the  ionosphere.  Therefore,  the 
multiple  scattering  process  must  be  considered  and  the  innate  nature  of  the 
quasi-particle  theory  is  particularly  suitable  for  the  study  of  ionospheric 
scintillation. 

The  wave  characteristics  including  its  amplitude  and  phase  can  be 
described  by  a  joint  distribution  function  of  the  position  vector  f  and  the 
wavevector  k ,  known  as  the  Wigner  Distribution  Function  which  is  defined 
to  be 

F(r,k,t)  =  7T^jdf'Qxp(2ik-r')i//*  (r  +r\t)y/{f -r\t)  (1) 

where  ^  is  a  wave  field  component  and  the  integration  is  over  all  space 
where  the  field  exists.  The  quasi-particle  spatial  density  and  spectral  density 
can  be  obtained  from  the  Wigner  distribution  function  by  the  two  relations 
|  ydr,t) |2  =  jF(T,k,t)dk  and  \G(k,t)\2  =  jF(r,k,t)df  ,  respectively. 

A  single-mode  wave  equation  is  considered  and  expressed  as 

i^jf  =  ~(y2  +ko2£\(r))V'  (2) 

The  corresponding  transport  equation  of  the  Wigner  distribution  function  is 
derived  [Ho  et  al.,  1994]  to  be  : 
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(3) 


k-VF  =  -(k02  /  ^)jdk'  F(r,k')jds{slr(f-s)sin[2(k-k'ys] 

+  £-ly(f-?)cos[2(* -A')  *]} 

where  £\{f)  =  sir)  -  £q  =  £\r(?)  +  isy(f)  and  e0  is  the  permittivity  of  the 
unperturbed  background. 

The  right  hand  side  of  the  transport  equation  (3)  can  be  written 

(k0+kx)-  +  ky—F  =  -(Ki/z1)Gi®F  (4) 

dx.  dy_ 

where 

G\{r,k  -  k')  =  J ds'{e^r{r  -  ;?)sin(2(£  -  k) •  s)+ey  (r  -  s)cos(2(£  -  k ')  *5)} 
and  <8>  stands  for  convolution. 

The  first  three  momentum  moments  of  the  quasi-particle  distribution 
function  are  p=jdkF(r,k),  J  =  jdkkF(r,k)  =^y?+Jj,  and 

s  =  j dkkkF(r,k)  =  %%?+  k$J\  +  J\kQ  +  ^ ;  which  are  the  quasi-particles' 

number  density,  current  density  and  stress  tensor  respectively.  The  governing 
moment  equations  can  be  derived  from  (4)  and  the  first  two  moment 
equations  are  found  to  be 

*0§+v-Ji=-vW^) 


li 


and 


(5) 


k°  dx *^1  +  ^  =1  _  ko2jleli (?)  +  ^o  2 ^e\r?) 

We  now  consider  wave  scattering  within  a  layer  of  plasma.  The 
plasma  is  perturbed  by  density  irregularities  represented  by  a  two  dimensional 
spatial  variation  of  the  form 

n{x,y)  =  Pp(x-  £){n0  +  Am[1  -  (sinfcj  jx  -sin&jQx)  /  Mjx} 

[\-(smk2iTj-smk2o7i)/ Ak2Tj]}  (6) 

where  rj  =  cos  6(y  +  x  tan  0) ,  Pf  {x)  =  1  for  \x\<£  and  0  for  |x|  >  £  represents 
a  gate  function.  This  density  variation  can  be  decomposed  into  two  separate 
components.  The  first  is  characterized  by  (ki)=(ku  +kio)/2=27t/di  with  the 
average  periodicity  dj  in  the  longitudinal  direction.  The  second  component 

has  an  average  periodicity  d2,  characterized  by  <k2>=(k21  +k2o)/2=27c/d2, 
oblique  to  the  first  with  an  angle  (/>=  n!2  -  0;  nQ  is  the  background  plasma 

density,  and  An  is  the  amplitude  of  the  density  irregularity. 

This  variation  consists  of  two  gratings  offset  from  each  other  by  an 
angle  ^  as  illustrated  in  Fig.  8.  As  shown  these  two  gratings  extend  to  infinity 

in  the  transverse  direction  and  have  a  finite  extent  along  the  downward 
propagation  direction. 

It  is  expected  that  a  model  provides  a  description  of  irregularities  and 
not  resulting  signal  scintillation  exclusively  [Fremouw  et  al.,  1984].  This 
model  provides  a  useful  description  of  BSS  irregularities  containing 
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Incoming  Plane  Wave 


Fig.  8.  Incoming  plane  wave  is  scattered  by  a  two-dimensional  irregularity 
grating  characterized  with  grating  wavevectors  (kj>  and  <k2>  offset 
by  angle  <t>.  The  thickness  of  the  irregularity  region  is  2 L 


information  on  irregularity  shape  as  well  as  strength.  In  addition 
parameterization  allows  adjustment  due  to  geophysical  changes. 

The  complex  dielectric  function  is  derived  from  the  plasma  density 
function  as 

e(x,y)  =  l -[&po*‘  v^)](l-/W c o)Pg(x-i ) 

=  s0  +  £l(x,y)  (7) 

where  Sj(g)  =  (smkn%-sinkj0%)/ }  kjQ  <  k.  <  k^  ?  Akj  =  k^  -  kiQ  for 

i=l,2.  S  is  the  Dirac  delta  function,  s0  =  1-  aPe(x-  t)(l-iv/  w), 
st  =  -Pa  Pe  (x  - 1) [1  -  Si (*)]  [1  -  &i ( rj)] ,  (Op  /  2  n  is  the  background  electron 

plasma  frequency,  CD  is  the  wave  angular  frequency,  v  is  the  electron-ion 
collision  frequency,  J3=An/tiQ  is  the  percentage  of  density  fluctuation  with 

respect  to  the  background  plasma  density,  and  a=  afp I (aP- +  v1)  « 1  is 
assumed. 

To  truncate  the  hierarchy  of  the  moment  equations,  a  cold  quasi- 
particle  distribution  is  assumed.  This  assumption  is  based  on  the  fact  that  the 
beacon  signal  is  monochromatic  and,  hence,  the  initial  quasi-particle 
distribution  has  zero  temperature.  Thus,  the  second  moment  tensor  =1 
becomes  zero  and  the  two  moment  equations  in  (5)  form  a  closed  set  of 
equations  for  further  analysis.  With  the  aid  of  (7)  the  two  equations  in  (5) 
become 
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k°  -£p+  V  •  Ji  =  a(3k^dr  v/  wP£  (x  -  f  )[1  -  8\  (*)][1  -  S2  (  7)]) 


and  (8) 

k°  & ^1 + v '  -1 =  aPio'~h  (~v/  wPe  (*  -  OP  -  (*)]P  -  si  ( ’/)]) 

+apk02p{  *f{[l  -  <J2 ( 7)]((<5(* - 20 -  <5(x))[l -  £,(*)] 

+  Pf,(x-  f  )[£j  (*)  /  JC  -  J  COSArj  I*  +  fcjQ  COS^jqX  /  MjX]) 
+P/(jf-^)[l~^(j:)](sm^/  rj)(82(rj) 

-^2i  cos^2jJc  +  &2q  cos^q^  /  A&2  ]) 

+5^(*-f)[l- <?i(x)](cos0/  rj)(82(?j) 

-^2i  cos^i*  +  ^20  cos^20*  ^  ^2  D 

The  Scintillation  Index  is  then  calculated  as: 

«S4=([</^2>-</^>2]  /  <P>2)  > 

where  the  average  is  performed  over  a  spatial  period  of  the  transverse 
irregularity,  only  the  intensity  distribution  p=\  y/(r,t)\2of  the  wave  is  needed. 

The  observations  showed  that  the  ionosphere  irregularities  having 
nearly  sinusoidal  waveforms  in  the  plasma  density  extended  up  to  7000  km  in 
the  E-W  directions.  The  rms  irregularity  amplitude  was  about  3.5%  of  the 
background  electron  density,  which  was  about  9.9x\$®m~^ .  The  thickness 
of  the  irregularity  layer  was  about  35  km.  The  satellite  was  located  at  an 
altitude  of  about  400  km  and  the  slant  range  of  the  signal  was  about  600  km. 
Thus,  in  the  numerical  analysis  we  have  chosen  Aki=0  and  dj=2^,  p=3.5%, 
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nQ=  9.9  x  10lOm-3,  2f=50km,  and  350km  for  the  free  space  propagation  from 
the  bottom  of  the  layer  to  the  ground. 

The  calculations  include  the  effect  of  the  finite  uniform  spectral  width 
Ak2/(k2)  of  the  density  irregularities  on  the  scintillation  index  S4  measured 
on  the  ground.  The  signal  frequencies  considered  are  250MHz  and 
1541MHz.  The  optimized  S4  is  evaluated,  for  each  frequency  by  choosing 
=(2^L)l/2  so  that  L=350km  is  the  Fresnel  size  of  the  signal  propagating 

down  from  the  exit  plane  of  the  irregularity  region  to  the  ground;  where  X  is 
the  signal  wavelength. 

Presented  in  Fig.  9  is  the  optimized  S4  for  250MHz  beacon  signal  with 
d2=  707  m,  an  additional  case  d2=  1000  m  is  also  shown.  The  dependence  of 
the  optimum  S4  on  spectral  width  Ak2/  (k2)  is  evaluated.  It  is  shown  that  the 
location  and  value  of  the  first  peak  of  the  S4  curve  does  not  varies  with  d2 
and  the  peaks  are  located  at  zero  spectral  width.  A  small  spectral  width  in  the 
irregularity  spectrum  reduces  the  optimized  S4  value,  the  tendency  of  the 
effect  of  spectral  width  is  to  reduce  the  S4  value.  It  shows  that  in  the 
optimized  case  S4  drops  to  a  minima  near  Ak2/(k2)  «  0.08.  As  Ak2/(k2) 
increases  further  above  0.08,  the  optimized  S4  curve  features  damped 
oscillation.  Since  the  values  of  S4  in  the  region  of  Ak2/(k2)  >  0.08  are  quite 
small,  Ak2/  <k2)<0.08  will  be  the  region  of  interest.  The  diminishing  effect  of 
the  spectral  width  on  the  S4  value  suggests  that  the  spectral  distribution  of 
the  irregularities  have  smeared  out  the  result  of  the  multiple  scattering  and 
hence  reduced  the  net  scintillation  effect  on  the  transverse  signals.  Similar 
functional  dependences  of  S4  for  a  1541MHz  signal,  corresponding  to  that 
used  by  the  geostationary  beacon  satellite  Fleetsatcom  [Basu  et  al.,  1987],  are 
presented  in  Fig.  10  showing  similar  features  results  of  Fig.  10  suggest 
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Ak2/(k2> 


Fig.  9.  Dependence  of  S4  on  spectral  width  Ak2/(k2)  for  250MHz  beacon 
signal,  with  <ij==50km,  707  m,  1000  m,  2f=50  km,  n0  =  1x10^ 

<j>=45°,  and  p  =  3.5%,  where  d2=707  m  corresponds  to  optimized  S4 


Ak2/(k2) 


Fig.  10.  Dependence  of  S4  on  spectral  width  Ak2/(k2)  for  1541MHz  beacon 
signal,  with  <ij=50km,  280  m,  500  m,  2f=50  km,  n0  —  lxiol  1 nr 3, 

(j)=45°,  and  p  =  3.5%,  where  d2=280  m  corresponds  to  optimized  S4 


that  S4  decreases  as  frequency  increases  which  is  consistent  with  common 
understanding.  In  addition,  it  is  shown  that  the  minima  location  of  S4  varies 
with  the  frequency  as  show  in  Figs.  9  &  10.  The  first  minimum  of  the 
optimized  S4  is  located  at  Ak2/(k2)=  0.08,  .02,  for  the  250MHz  and 
1541MHz  cases  respectfully.  The  separation  of  the  locations  of  the  S4 
minima  for  different  beacon  signal  frequencies  could  be  useful  for 
distinguishing  the  mechanism,  (  irregularity  amplitude  or  large 
spectralwidth),  responsible  for  small  S4  values. 

The  results  of  Fig.  10  suggest  that  S4  decreases  as  frequency 
increases  which  is  consistent  with  common  understanding.  In  addition,  it  is 
shown  that  the  minima  location  of  S4  varies  with  the  frequency  as  show  in 
Figs.  9  &  10.  The  first  minimum  of  the  optimized  S4  is  located  at  Ak2/(la)= 
0.08,  .02,  for  the  250MHz  and  1541MHz  cases  respectfully.  The  separation 
of  the  locations  of  the  S4  minima  for  different  beacon  signal  frequencies 
could  be  useful  for  distinguishing  the  mechanism,  ( irregularity  amplitude 
or  large  spectral  width),  responsible  for  small  S4  values. 
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ABSTRACT 

In  the  magnetosphere  energetic  charged  electrons  in  the  radiation  belts  are  trapped 
by  Earth's  magnetic  field.  In  the  equatorial  region  where  a  symmetric  mirror  field  may  be 
assumed,  the  electrons  undergo  bounce  motion  along  the  lines  of  force.  When  a  large 
amplitude  whistler  wave  is  present,  the  motion  of  the  electrons  becomes  perturbed.  The 
nonlinear  interaction  due  to  the  spatial  dependence  of  the  field  quantities  causes  the 
motion  of  some  of  the  trapped  particles  to  become  chaotic.  This  chaotic  behavior  serves 
to  enhance  the  electrons'  axial  kinetic  energy,  scattering  some  electrons  into  the  magnetic 
mirror's  loss  cone.  Thus  we  offer  a  possible  explanation  of  the  observed  electron 
precipitation  in  the  polar  regions. 

Introducing  a  canonical  transformation,  the  Hamiltonian  of  a  trapped  electron 
becomes  time  independent.  With  the  aid  of  additional  constants  of  motion,  we  reduce  the 
equations  of  motion  so  that  the  electron  trajectories  lie  in  a  three  dimensional  phase  space. 
We  then  use  a  Surface  of  Section  technique  to  examine  the  chaoticity  of  the  system. 
Lyapanov  exponent  analysis  is  also  used  to  extend  our  results  to  the  case  of  very  small 
bounce  frequency,  as  one  would  expect  in  the  magnetosphere.  Pitch  angle  scattering 
techniques  are  also  effective  in  describing  the  system  with  small  bounce  frequency,  and 
confirm  the  enhancement  of  the  electron  axial  kinetic  energy. 


INTRODUCTION 


In  the  magnetosphere,  charged  particles  (electrons  and  hydrogen  ions)  are  trapped 
by  the  Earth's  magnetic  field.  When  these  particles  are  undisturbed  they  undergo  bounce 
motion  about  the  geomagnetic  equator.  However,  if  a  large  amplitude  plasma  wave  is 
present  in  the  magnetosphere,  particle  trajectories  may  become  significantly  perturbed. 
This  allows  some  particles  to  be  scattered  into  the  mirror  loss  cone,  and  precipitate  into 
the  polar  regions  along  the  geomagnetic  field  lines  where  they  have  been  observed 
experimentally. 

In  the  case  of  electrons,  the  original  experimental  evidence  for  the  precipitation 
was  in  the  observation  of  "Trimpi  events",  named  for  their  discoverer  [Helliwell,  et.  al., 
1973],  These  events  are  the  perturbations  of  amplitude  and/or  phase  of  VLF  and  LF 
beacon  signals  which  are  trapped  in  the  Earth  ionosphere  waveguide  and  are  used 
primarily  for  navigation  purposes.  An  event  is  characterized  by  a  rapid  onset  (~1  s)  of  the 
perturbation  in  the  detected  signal  and  a  slow  recovery  (-100  s)  back  to  the  original 
conditions.  These  perturbations  are  related  to  secondary  ionization  in  the  ionospheric  D 
region  caused  by  energetic  electrons  precipitating  from  the  magnetosphere  [Chang  and 
Inan,  1983].  These  secondary  ionizations  alter  the  local  plasma  density,  and  thus  the 
dielectric  constant  of  the  Earth  ionosphere  waveguide.  This,  in  turn,  accounts  for  the 
variation  in  the  beacon  signal's  received  amplitude  or  phase,  as  the  onset  of  these 
precipitation  events  is  rapid  but  the  return  to  the  quiescent  state  is  determined  by  the 
longer  atmospheric  recovery  time.  More  recent  experiments  have  given  direct 
confirmation  of  the  presence  of  energetic  electron  precipitation  [Amoldy  and  Kinter, 
1989],  In  these  experiments,  electron  spectrometers  have  been  carried  aloft  either  in 
satellites,  high  altitude  balloons,  or  atop  rockets. 

Common  to  all  methods  of  observing  the  electron  precipitation  is  the  one  to  one 
correspondence  of  the  presence  of  magnetospheric  whistler  waves  with  the  precipitation 
event  [Chang  and  Inan,  1985].  Whistler  waves  have  been  observed  propagating  between 
the  Northern  and  Southern  hemispheres  along  magnetic  field  lines.  These  waves  can  be 
coupled  into  the  magnetosphere  either  from  natural  sources,  such  as  lightning  [Chang, 
1985],  or  man  made  ones,  such  as  high  power  VLF  waves  transmitted  for  ionosphere 
modification  experiments  [Wright,  1975].  In  the  case  of  lightning,  the  sudden  current 
impulse  produces  radiation  throughout  the  VLF  range.  These  waves  are  normally  trapped 
by  the  Earth  ionosphere  waveguide.  However,  the  presence  of  numerous  "ducts",  or  field 
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aligned  ionization  density  irregularities  can  trap  electromagnetic  radiation  and  guide  it 
between  hemispheres  [Carpenter  and  Sulic,  1988].  The  existence  of  these  ducts  is 
confirmed  by  the  ground  based  observation  of  whistlers  which  do  not  undergo  total 
reflection  upon  reaching  the  ionosphere-magnetosphere  boundary  [Helltwell,  1965]. 

Due  to  the  correlation  of  the  presence  of  magnetospheric  whistler  waves  and  the 
electron  precipitation  events,  there  has  been  much  theoretical  work  on  the  interaction  of 
the  waves  and  the  trapped  radiation  belt  electrons.  All  of  the  previous  work  has  focused 
on  the  interaction  of  waves  whose  frequency  matches  that  of  the  electrons'  Doppler  shifted 
cyclotron  resonance  frequency  [Kennel  and  Petschek,  1966,  Inan,  1987,  Villaldn  and 
Burke,  1991],  a>-k  u  =  Q.,/y,  where  y  =  (l-t?/c2 )'1/2  is  the  relativistic  factor.  While 
this  theory  has  been  successful  in  explaining  some  of  the  events,  there  is  experimental 
evidence  that  other  processes  are  also  important  [Goldberg,  et.  al.,  1987]. 

Additional  evidence  that  this  theory  may  require  augmentation  derives  from  the 
simultaneous  observation  of  precipitation  events  at  geomagnetically  conjugate  regions  in 
both  hemispheres  due  to  a  single  lightning  flash  [Burgess  and  Inan,  1990].  In  order  to 
match  the  resonance  condition,  particles  must  move  in  the  opposite  direction  as  the  wave, 
since  the  wave  frequency  co  is  smaller  than  the  electron  cyclotron  frequency  Cl,  for  the 
whistler  wave.  Thus  this  process  would  first  precipitate  electrons  back  to  the  launching 
point  of  the  wave.  In  order  induce  precipitation  into  the  opposite  hemisphere,  either  the 
electrons  themselves  must  be  reflected  by  the  magnetosphere-ionosphere  boundary  and 
bounce  back  to  the  opposite  hemisphere,  or  the  whistler  .wave  must  be  reflected  at  the 
boundary  and  interact  with  particles  moving  in  the  other  direction  on  its  return  path.  In 
either  event,  there  would  be  an  appreciable  delay  between  precipitation  events  in  opposite 
hemispheres,  due  to  the  time  required  for  the  particles  or  wave  to  bounce  from  one 
hemisphere  to  the  other.  Thus  we  would  like  to  propose  an  alternate  theory  for  the 
electron  precipitation  that  does  not  favor  one  direction  of  precipitation  over  the  other,  and 
therefore  is  consistent  with  the  simultaneous  observation  of  precipitation  events  in  both 
hemispheres.  Shown  in  Fig.  1  is  an  overview  of  the  situation  considered. 

Previous  work  has  studied  the  nonlinear  interaction  between  the  magnetospheric 
protons  and  the  axial  component  of  the  kinetic  Alfven  wave  (KAW)  fields  with 
applications  to  particle  precipitation  [Prakash,  1989],  It  was  found  that  a  large  amplitude 
KAW  will  induce  chaos  in  a  test  particle's  orbit .  Further  studies  [Ho,  et.  al.,  1994]  found 
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that  the  inclusion  of  the  transverse  components  of  the  wave  fields  not  only  lowers  the 
wave  amplitude  threshold  for  the  onset  of  chaos,  but  the  chaos  also  enhances  the  particle's 
transverse  (radial)  displacement,  allowing  the  precipitation  of  some  particles  into  the 
equatorial  region  [Guzik,  et.  al.,  1989],  in  addition  to  the  loss  cone.  Confirming  evidence 
for  these  theories  has  also  been  provided  by  a  recent  chamber  experiment  [Warren  and 
Mauel,  1995],  Since  chaos  is  believed  to  be  important  to  the  magnetospheric  phenomena, 
and  whistler  wave-electron  interactions  cause  the  precipitation  events,  we  continue  our 
work  [Faith,  et.  al.,  1995]  on  the  dynamics  of  a  single  trapped  electron  subjected  to  a 
large  amplitude  whistler  wave. 


FORMULATION 

To  simplify  the  formulation  while  retaining  all  the  essential  physics  the  magnetic 
dipole  field  is  replaced  by  a  parabolic  scalar  potential,  cp  oc  z2,  superimposed  over  a 
uniform  magnetic  field,  B  =  .  The  parabolic  potential,  characterized  by  a  bounce 

frequency  co6,  simulates  the  mirror  effect  of  the  dipole  field.  Justification  for  this 
simplification  is  as  follows: 


A  parabolic  mirror  field  in  cylindrical  coordinates  may  be  describes  as: 

B  =  fr0[l  +  (z2  -r2)  /  L2]z  -  BqZ?  /  L2 , 

where  L  is  the  scale  length  of  the  magnetic  field.  One  can  show  that  this  magnetic  field 
profile  satisfies  both  V*  Z?  =  0  and  V  xB  =  0,  as  one  would  expect  in  the  magnetosphere. 
The  z  component  of  the  equation  of  motion  of  an  electron  in  the  above  field  is  given  by: 


— v  =  —z  =  ~—Bnzrv  11}  =  -- E3L*- 
dt  '  dt2  m  ^ 6 mL 


where  p  =  ervQ  /  2  is  the  magnetic  dipole  of  the  electron,  and  a>b  =  ([iB(/nL2)V2  is  the 


bounce  frequency.  Thus  a  parabolic  potential  may  be  used  to  simulate  the  mirroring  effect 
of  the  magnetic  dipole  field.  The  simplification  cofc  =  const,  is  proper  as  long  as  the 

magnetic  moment  varies  slowly  in  time.  It  is  also  noted  that  the  approximation  B  =  B£ 
applies  only  when  z2  /  Z?  « 1  and  r2  / 1}  « 1 . 


The  whistler  wave  field  may  be  described  by  the  vector  potential: 

\  =  (.5 /^[fcos(fc-©0+.P  sin  (&-©*)]•  (1) 
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The  total  vector  potential  has  contributions  from  both  the  wave  and  static  fields, 
and  is  given  by  A  =  Aw+yB0x.  The  required  scalar  potential  is  given  by 
<p  =  -(1  /  2e)m®2bz2 .  With  both  of  the  potentials  given  we  can  derive  the  classical  electron 
Hamiltonian  as: 

H  =  (\/2m)(p+eAf  -eq> 


+^Lj[/?x  cos (Az  —to/) -mOx sin(Az  -©/)]  +/?^£2x-^Lsin(fc  -©/)j  (2) 

where  Q  =  eB0/ m,  and  -eBIm  are  the  electron  cyclotron  frequency  and  a  frequency 
proportional  to  the  wave  magnetic  field  amplitude,  respectively. 


The  equations  of  motion  corresponding  to  the  above  Hamiltonian  will  give 
trajectories  in  a  six  dimensional  (plus  time)  phase  space.  With  a  phase  space  of  this  many 
dimensions  it  is  difficult  to  present  results  that  have  clear  physical  meaning.  We  therefore 
seek  to  simplify  the  Hamiltonian  in  a  number  of  ways.  First,  since  energy  (equal  to  the 
Hamiltonian)  is  defined  only  to  within  an  additive  constant,  the  constant  term 
{ml2){D^I  kf  may  be  neglected  without  altering  the  system  dynamics,  and  is  dropped. 
Next,  the.y  coordinate  is  found  to  be  cyclic,  meaning  py  is  a  constant  of  motion  which  can 

be  set  to  zero  without  loss  of  generality.  Finally  we  wish  to  eliminate  the  time  dependence 
from  the  Hamiltonian. 


We  accomplish  this  by  introducing  the  generating  function. 


=  \m£l{x2  +  (T2)  cot©/ -mflxQ  esc©/ . 


(3) 


This  function  provides  a  canonical  transformation  with  new  canonically  conjugate 
coordinates.  The  new  coordinates  (Q,P)  are  related  to  the  old  (x,px)  through  the 
generating  function  by: 


which  gives: 


P  =  - 


a 5. 

dQ’ 
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P  =  px  cos((i)t)+mQxsm((i)t) 

Q  =  -(px/mQ)sm(d)t)  +  xcos((s>t). 


(4) 


In  addition,  we  can  derive  the  new  time  independent  Hamiltonian  from  the  old  one  and  the 
generating  function  through: 


dt 


+— (Pcosfe-wQQsinAz)  (5) 

k 

Using  Hamilton's  equations  of  motion: 
r  =  VpK,  p  =  -VAT, 

and  the  normalizations  to  dimensionless  quantities: 


kQ  ->  Q,  kP  / /n©  — »  P,  kz  -»  z,  kpjm co  — >  pz, 
n/©->n,  and  ©6/©-^©fc, 


we  can  derive  the  final  set  of  autonomous  equations  of  motion  to  be  integrated 
numerically: 

£)  =  [(Q  - 1)  /  ft]/1 + J2,  cosr 
P  =  -D(Q-  1)^+Qn,  sin  z 
i  =  P: 

pz  =  -©jr+QI(/,sinz+flgcosz)  (6) 


Since  K  is  a  constant  of  motion,  the  set  of  equations  (6)  describes  trajectories  in  a 
three  dimensional  phase  space.  Thus  a  surface  of  section  technique  can  be  used  to  reduce 
the  three  dimensional  continuous  time  system  to  a  two  dimensional  map,  where  we  can 
examine  the  chaoticity  of  the  system  graphically.  It  should  be  stressed  that  the  bouncing 
motion  introduces  the  chaotic  nature  to  the  system.  One  can  show  that  if  there  is  no 
bouncing  motion,  i.e.  ©6  =  0,  that  an  additional  constant  of  motion: 

j(n2Q2  +  P2)-Opz  =  const. 
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can  be  found.  This  is  easily  verified  with  the  use  of  the  Poisson  Bracket.  This  constant 
reduces  the  system  to  two  autonomous  equations.  A  system  of  this  type  can  never  exhibit 
chaotic  motion,  and  the  trajectory  will  always  have  regular  behavior  in  a  two  dimensional 
phase  space. 


NUMERICAL  ANALYSIS 

In  the  following  analysis,  with  cob  *  0,  we  choose  to  eliminate  the  variable  P  with 
the  constant  of  motion  K.  We  choose  Q  =  0  as  the  surface  of  section,  allowing  us  to 
graphically  depict  the  chaoticity  of  the  system  with  a  portrait  of  the  (z,p,)  phase  space.  It 
is  noted  that  for  a  given  K,  both  P  and  Q  are  double  valued.  Hence  the  surface  of  section 
is  separated  into  the  Q<  0  and  Q>  0  cases.  However,  since  the  two  cases  are  mirror 
images  of  each  other  we  choose  only  to  present  the  Q  <  0  case. 

The  plots  are  made  by  integrating  the  equations  of  motion  forward  in  time  using  a 
fifth  order  Runge  Kutta  ordinary  differential  equation  solver.  When  the  sign  of  Q  changes 
from  positive  to  negative  between  two  time  steps,  the  trajectory  between  these  two  points 
is  interpolated  to  the  Q  =  0  plane,  and  the  resulting  point  recorded.  The  time  integration 
continues  for  many  (-1000)  bounce  periods.  The  several  orbits  on  each  of  the  surface  of 
section  plots  are  computed  by  beginning  multiple  time  integrations  from  different  initial 
conditions.  These  are  chosen  by  maintaining  a  constant  transformed  Hamiltonian  K,  and 
choosing  a  sequence  of  P=nA,  where  n=0,l,2,...  is  an  integer  and  A  is  a  value  chosen  to 
give  about  twelve  sets  of  initial  conditions.  These  values  of  P  are  used  to  solve  for  the 
proper  pz  (as  long  as  pz  remains  real!).  We  find  that  the  behavior  of  the  system  is  strongly 
dependent  on  two  quantities,  the  value  of  K  and  the  wave  amplitude.  As  flj  is 
proportional  to  the  wave  amplitude  and  K  equals  the  initial  particle  energy,  we  vary  these 
parameters  from  small  to  relatively  large  values  to  demonstrate  the  transition  from  regular 
behavior  to  chaos. 

For  our  analysis,  we  choose  parameters  corresponding  to  conditions  in  the 
magnetosphere,  with  the  exception  of  the  bounce  frequency  which  is  initially  chosen  much 
larger  than  the  physical  bounce  frequency  to  allow  the  use  of  more  modest  computation 
times  and  data  file  sizes.  The  results  will  then  be  extended  to  the  range  of  bounce 
frequencies  one  would  physically  expect.  The  other  parameters  are  chosen  to  reflect  a 
typical  whistler  wave  propagating  along  the  L=3  shell,  where  QJ2k=32  kHz,  co/27t=13.7 
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kHz,  and  a  plasma  frequency,  a>  p  1 2n  =  (ne2  /  zjri)m -15kHz.  The  plasma  frequency  is 
required  to  find  the  wavenumber,  k,  which  is  related  to  ©  through  the  whistler  wave 
dispersion  relation  ©  =  Cl,c2k2  /©2.  As  a  starting  point  we  also  choose  the  bounce  and 

wave  frequencies  equal. 

First  we  choose  to  fix  the  particle  energy  and  vary  the  wave  amplitude.  K  is 
chosen  to  be  1,  which  corresponds  to  the  commonly  observed  particle  energy  of  180  keV. 
With  a  very  small  wave  amplitude,  Q,  =  0.1,  the  trajectories  in  Fig.  2a  are  found  to  be 
regular  for  all  initial  conditions.  Increasing  the  wave  field  to  fl,  =  0.4  in  Fig.  2b,  the 
trajectories  are  noticeable  distorted  by  the  larger  wave,  and  the  entire  surface  of  section  is 
pushed  into  the  pt<  0  portion  of  phase  space.  Continuing  to  increase  the  wave  field  to 
Qj  =  0.43  for  Fig.  2c  the  trajectories  for  several  initial  conditions  become  chaotic,  as 
shown  by  the  irregular  smear  of  dots  in  the  upper  portion  of  the  plot.  Finally,  one  last 
increase  in  the  wave  amplitude  to  Q,  =  0.6  in  Fig.  2d  shows  all  of  the  regular  orbits 

replaced  by  chaos  throughout  the  plot. 

The  second  parameter  that  strongly  influences  the  dynamics  of  the  system  is  the 
particle  energy.  For  all  the  plot  in  this  figure,  the  cyclotron  ,  wave,  and  bounce 
frequencies  are  all  held  constant  at  Q  =  2.333,  Q,  =  0.46,  and  ©t  =  1.0.  In  the  first  plot, 
Fig.  3a,  the  initial  particle  energy  is  45keV,  corresponding  to  K=0.25.  The  trajectories  are 
regular  for  all  initial  conditions.  With  K=0.5,  Fig.  3b,  the  trajectories  all  remain  regular, 
although  some  distortion  is  seen  from  the  previous  case.  Once  the  particle  energy  is 
increases  to  K=0.7,  Fig.  3c,  the  onset  of  chaos  has  been  observed  in  the  heavy  smear  of 
dots  in  the  center  of  the  plot.  Finally,  a  small  further  increase  in  K  to  K=0.75,  Fig.  4d, 
caused  the  trajectories  to  become  chaotic  over  the  majority  of  the  plot. 

To  extend  the  results  to  a  more  realistic  bounce  frequency,  10'3  £©*  £  10"2,  two 
techniques  are  available  to  us.  The  first  involves  the  use  of  Lyapanov  exponent  analysis. 
These  exponents  are  the  rates  of  exponential  separation  of  two  initially  close  points  in 
phase  space.  As  this  exponential  separation  defines  a  chaotic  system,  a  positive  Lyapanov 
exponent  means  the  system  is  chaotic.  Presented  in  Fig.  4  is  the  dependence  of  the  first 
two  Lyapanov  exponents  on  the  bounce  frequency.  Since  in  Hamiltonian  systems  the 
exponents  occur  only  in  oppositely  signed  pairs,  only  the  positive  exponents  are  presented. 
The  exponents  were  computed  for  five  bounce  periods  using  the  method  of  Wolf,  et.  al., 
[1985],  This  time  is  sufficient  for  the  exponents  to  settle  to  their  asymptotic  values. 
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The  second  way  to  extend  our  results  to  the  case  of  small  bounce  frequencies  is  to 
define  the  pitch  angle  a  as: 

a  =  tan-'  (v±  /  v,)  =  tan-1  (^/vj+vj  /  vx),  (7) 

where  the  true  velocities  are  related  to  the  canonical  momenta  and  our  transformed  and 
normalized  variables  as: 

mvx  =  px+eAx  =  P  cos©/  -  mQQ  sin  ©/  +  (mQi  /  k)  cos (©/ -kz) 

— »  P  cost  -  CIQ  sin  / + cos(/  -  z) 

mvy  =  Py+eAy  =  Psin©/+mQ0cos©/  +  (/nQ1  /£)sin(©/-Az) 

->  Psin/+Q0cos/  +  n,  sin(/-z) 


mvz  =  Pz~>P:- 

The  quantity  a  is  especially  of  interest  when  the  electron  is  at  the  geomagnetic  equator 
(z=0),  where  it  can  be  compared  with  the  loss  cone  angle.  If  at  this  plane  the  pitch  angle 
is  less  than  the  loss  cone  angle  the  particle  has  enough  axial  momentum  to  escape  out  the 
end  of  the  magnetic  mirror  field,  and  precipitate  into  the  polar  regions.  Shown  in  Fig.  5  is 
the  equatorial  pitch  angle  plotted  against  time  for  four  bounce  frequencies, 
ab  =  1.0,10“',10"2,10"3,  in  Figs.  5a,  5b,  5c,  and  5d,  respectively.  In  all  cases  it  is  found 
that  in  the  chaotic  case  the  pitch  angle  varies  wildly  and  for  all  bounce  frequencies 
achieves  lower  pitch  angles,  and  can  more  readily  escape. 

CONCLUSION 

We  have  shown  that  the  nonlinear  interaction  arising  from  the  spatial  dependence 
of  the  fields  of  a  large  amplitude  whistler  wave  with  a  trapped  bouncing  electron  can  lead 
to  chaos  in  the  electron  trajectory.  In  addition,  we  have  found  that  it  is  the  bouncing 
effect  that  leads  to  the  chaos,  as  shown  by  the  existence  of  an  additional  constant  of 
motion  in  the  case  of  zero  bounce  frequency.  However,  it  should  be  noted  that  our  model 
of  the  geomagnetic  field,  that  is  of  a  parabolic  potential  superimposed  over  a  uniform 
magnetic  field,  is  valid  only  in  the  region  near  the  geomagnetic  equator. 
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The  end  result  of  the  chaos  induced  in  the  electron  trajectory  is  to  enhance  its  axial 
kinetic  energy,  as  shown  in  any  of  the  chaotic  surface  of  section  plots  in  the  (z,pt)  phase 
space.  This  is  apparent  because  when  the  electron  trajectory  is  regular,  it  is  confined  to 
lay  only  upon  the  closed  loops  seen  in  Fig.  2a  or  Fig.  3a.  When  the  wave  amplitude 
becomes  larger,  and  the  trajectories  become  chaotic,  the  electron  trajectory  is  no  longer 
confined  to  a  closed  loop,  and  may  wander  throughout  the  chaotic  regions.  Thus  an 
electron  with  an  initially  small  axial  kinetic  energy  may  eventually  acquire  enough  energy 
to  be  scattered  into  the  loss  cone  of  the  magnetic  mirror  field.  In  addition,  this  increase  in 
axial  kinetic  energy  will  alter  the  pitch  angles  of  the  trapped  particles,  as  seen  in  Fig.  5., 
and  allow  scattering  into  the  loss  cone. 

Once  electrons  are  scattered  into  the  loss  cone,  they  follow  the  geomagnetic  field 
lines  and  will  precipitate  into  the  ionosphere,  where  they  can  be  detected  by  the  methods 
discussed  in  the  introduction  (and  in  the  references  therein).  The  dissimilarity  of  this 
theory  as  compared  to  one  based  on  gyroresonant  interaction  is  the  fact  that  here  there  is 
no  preferred  direction  of  precipitation.  Thus  particles  are  equally  likely  to  be  scattered  in 
the  direction  with  or  against  the  wave,  it  offers  a  possible  explanation  of  the  simultaneous 
observation  of  precipitation  events  in  geomagnetically  conjugate  regions  in  both  the 
Southern  and  Northern  hemispheres  [Burgess,  1990],  This  leads  us  to  believe  that  the 
chaotic  scattering  discussed  here  can  successfully  augment  the  gyroresonant  scattering 
theories. 
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Figure  2  (c)  and  (d) 

Surface  of  section  plots  at  the  Q=0  plane,  with  K=l,  H  =  2.333,  ob  =  1.0,  and  (c) 
Qj  =  0.43  showing  both  chaotic  and  regular  regions,  and  (d)  O,  =  0.6  showing  only 
chaotic  regions. 


Surface  of  section  plots  at  the  Q= 0  plane,  with  Q  =  2.333,  Q,  =  0.46,  and 
ob  =  1.0,  and  (a)  K=0.25,  and  (b)  K=0.5.  Both  plots  show  only  regular  behavior. 


(d) 

Figure  3  (c)  and  (d) 

Surface  of  section  plots  at  the  Q= 0  plane,  with  Q  =  2.333,  fi^O.46,  and 
(ob  =  1.0,  and  (c)  K=0.7,  and  (d)  K=0.75.  Both  plots  show  both  regular  and  chaotic 
behavior. 


Lyapanov  Exponents 


Figure  4  ^ 

Functional  dependence  of  first  two  Lyapunov  exponents  an  the  bounce  frequency, 

K=l,  Q  =  2.333,  =  0.46,  &  =  *o  =  0,  P0  =  l.2,  Ao  = °-27>  311(1  10“3  <cofc  <  i.  The 

system  is  chaotic  for  all  bounce  frequencies. 


Pitch  Angle  (Degrees)  Pitch  Angle  (Degrees) 


Figure  5  (a)  and  (b)  .  -  v_. 

Functional  dependence  pitch  angle  at  the  geomagnetic  equator  on  time  tor  K  l, 

Q  =  2.333,  Qo~z o  =  0>  =  1.0  (a),  and  cofc  =  0.1  (b).  The  chaotic  plot  is  for 

nj  =  0.46,  Pq  =  1.6,  and  /?z0  =  0.46,  while  the  regular  is  for  ft,  =  0.1,  P0  =  12,  and 

pz0  =  0.27 .  Both  plots  have  a0  *  70° . 
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Figure  5  (c)  and  (d) 

Functional  dependence  pitch  angle  at  the  geomagnetic  equator  on  time  for  K=l, 
ft  =  2.333,  Qo  =  z0  =  0,  d)b  =  \0~2  (c),  and  ©b  =  10-3  (d).  The  chaotic  plot  is  for 
Qj  =  0.46,  P0  =  1.6,  and  pz0  =  0.46,  while  the  regular  is  for  Q,  =  0.1,  P0  =  1.2,  and 

pz0  =  0.27 .  Both  plots  have  a0  «  70°. 
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